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ABSTRACT 

A star around a massive black hole can be disrupted tidally by the gravity of the 
black hole. Then, its debris may form a precessing stream which may even collide with 
itself. In order to understand the dynamical effects of the stream-stream collision on 
the eventual accretion of the stellar debris onto the black hole, we have studied how 
gas flow behaves when the outgoing stream collides supersonically with the incoming 
stream. We have investigated the problem analytically with one- dimensional plane- 
parallel streams and numerically with more realistic three-dimensional streams. A 
shock formed around the contact surface converts the bulk of the orbital streaming 
kinetic energy into thermal energy. In three-dimensional simulations, the accumulated 
hot post-shock gas then expands adiabatically and drives another shock into the low 
density ambient region. Through this expansion, thermal energy is converted back to 
the kinetic energy associated with the expanding motion. Thus, in the end, only a 
small fraction of the orbital kinetic energy is actually converted to the thermal energy, 
while most of it is transferred to the kinetic energy of the expanding gas. Nevertheless 
the collision is effective in circularizing the debris orbit, because the shock efficiently 
transforms the ordered motion of the streams into the expanding motion in directions 
perpendicular to the streams. The circularization efficiency decreases, if two colliding 
streams have a large ratio of cross sections and a large density contrast. But even in 
such cases, the main shock extends beyond the overlapping contact surface and the 
high pressure region behind the shock keeps the stream of the larger cross section from 
passing freely. Thus the stream-stream collisions are still expected to circularize the 
stellar debris rather efficiently, unless the ratio of the cross sections is very large (i.e., 
(Tl/(T 2 > 10). 

Subject headings: hydrodynamics - shock waves - galaxies : nuclei 
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1. Introduction 

It is believed that the tidal disruption of stars can occur at a rate of approximately one in 
10 4 years in typical nuclei of nearby galaxies with observational evidence for a super massive black 
hole (e.g., Rees 1988; Cannizzo, Lee & Goodman 1990). A large amount of energy (~ 0.1 m*c 2 ) 
can be released if a significant fraction of stellar mass is accreted onto the black hole. However, 
galaxies with convincing dynamical evidence for a central black hole usually do not reveal strong 
activities in their central parts (i.e., large M/L). It is possible that the dynamical mass implied 
by spectroscopic observations might not be due to the massive black hole but a cluster of dark 
stars (Goodman & Lee 1989). High resolution observations in the near future would enable us to 
eliminate some candidates for the constituents of dark clusters. It is, however, too early to rule 
out the existence of a central massive black hole on the grounds of low-level activities in galactic 
nuclei. We simply might not understand well enough the physical processes of the tidal disruption 
of stars and the subsequent evolution of the debris. 

One of the key questions here is how the stellar debris forms an accretion disk around the black 
hole and is accreted onto it (Rees 1994). Since the initial orbit of the debris is extremely eccentric, 
the accretion time scale is longer than the typical interval between two successive tidal disruption 
events (Gurzadyan & Ozernoy 1980, 1981). But, if the orbit of the debris can be circularized, the 
accretion time scale can become as short as a few tens of years (Cannizzo, Lee & Goodman 1990). 
It has been shown that strong relativistic effects make the orbit of the debris stream precess at a 
large rate around the black hole (on the order of one radian per revolution: see, for example, Rees 
1988; Cannizzo, Lee & Goodman 1990; Monaghan & Lee 1994; Kochanek 1994). As a result, the 
parts of the stream with different orbital periods (we will call them the outgoing and incoming 
streams, hereafter) may cross each other, and this stream-stream collision may dissipate the orbital 
energy via a shock leading to an effective circularization. 

So the circularization process through the stream-stream collision is a crucial issue in estimating 
the observational consequences of the tidal disruption of a star. However, the supersonic collision 
of two gas streams has not been investigated in detail, in order to tell us how the gas flow evolves 
subsequently and how accretion follows. One can expect that the details of the stream-stream 
collision and the efficiency of the circularization will depend upon the geometry and flow parameters 
of the gas streams at the time of the collision. Kochanek (1994) estimated that the widths and 
heights of the two colliding streams would be different by a factor of two or more, and the collision 
angle between the streams is typically 130° to 140°. He suggested that most of the larger stream 
would freely move away along its original path, because of significantly different cross sections, so 
that the stream-stream collision should have only small effects on the orbital evolution. 

Motivated by this, we have studied the evolution of the collision of the outgoing and incoming 
streams and estimated the degree of the conversion of the stream's orbital kinetic energy into other 
forms of energy. First, we have studied analytically the collisions of two one-dimensional, semi- 
infinite streams and showed how the shock parameters depend on the stream density and velocity. 
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In more realistic situations, the cross sections of the two colliding streams can be very different and 
also the collision is not necessarily head-on. Thus we have performed three-dimensional numerical 
simulations for the collisions of two supersonic streams with different density ratios, cross section 
ratios, and collision angles. 

In §2 we present the one- dimensional analytic study. This is followed by the description of 
the three-dimensional numerical study in §3. In §4 we describe the results and implications of the 
simulations. Conclusions are given in §5. 



2. One-dimensional Plane-Parallel Collision 

In this section, we describe the evolution of the collision of two plane-parallel gas streams, for 
which we can find a simple solution. It tells us some essential features found in the three-dimensional 
simulations. 

Suppose that two one-dimensional, semi-infinite, gas streams have density, p\ and p 2 , pressure, 
Pi and p 2 , and velocity, v i > and v 2 < and that they collide at / = 0. This is an example of the 
Riemann problem where a gas flow containing an arbitrary discontinuity sets into motion. Fig. 1 
shows the solution of the case where p\ = 4, p 2 = 1, p 1 = p 2 = 2.4 X 10~ 6 , and v\ = — v 2 = 0.2. 
The adiabatic index is 7 = 5/3. The solution for the post-shock regions 3 and 4 can be found from 
the intersection of the two Hugoniot curves on the so-called (p, v) diagrams for the two shocks (see, 
e.g., Zel'dovich & Raizer 1966). Here, the two shocks propagate from both sides of the contact 
discontinuity (CD). The pressure and velocity are continuous across the contact discontinuity, i.e., 
p 3 = p 4 and V3 = V4 = vcDi where vqd 1S the velocity of the contact discontinuity. 

Since the shocks are strong (i.e., p 3 , p^ ^ p\, P2, so P3/ pi = Pa/ P2 = (7 + l)/(7 — we nn( i 
the following simple relations: 



7 + 1 P1P2 , , 2 
P3 = Pa = , r- , — ^=^( v i ~ v 2), 



V 3 = V 4 = V C D 



Pi + VP2 



7 + 1 yp2 i . 

2 Jpi + Jp2 



7 + 1 ypi ( s. m 
° 2 = — TWVjw^ - " 2) (4) 



Here u\ and u 2 are the shock velocities relative to the upstream flows in the regions 1 and 2, 
respectively. The ratio of the shock velocities is |wi/w2| = \J P2/P1, so the ram pressures by the two 
shocks are same (i.e., p\u\ = p 2 u 2 ). 
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One can define the thermal energy conversion flux, Fth, as the postshock thermal energy 
generated by the shock per unit area per unit time according to 

F th = J^T^( W 2 - wi). (5) 

Here w\ = u\ + v\ and w 2 = u 2 + v 2 are the shock velocities in the frame where the initial flows 
are defined. In the strong shock limit, 

J? 7 + 1 PlP2 ( S 3 /RS. 



Then, the conversion ratio, R, can be defined by the ratio of the thermal energy conversion flux to 
the influx of the kinetic energy; that is, 

Fth 

R= (l/2) Pl vf + (l/2)p 2 vl (7) 

The probable configuration for the collisions of the debris streams mentioned in §1 is the one 
with x = /°i//°2^1 and v 2 ~ — v i- In this case, |mi/m2| = y/x- Then the postshock pressure is given 

by 

7+1 , s2 
P3 = -^—P2{Vl ~ V 2 ) 



(v^+i) 2 , 

if the shocks are strong. The term x/(\/x + I) 2 varies from 1/4 to 1 for x ^ 1- Thus, the postshock 
pressure is determined mostly by the relative velocity of the gas streams and the density of the 
lower density stream, but it is almost independent of the density of the higher density stream. 
Since F t h oc p^{vi — v 2 ), the thermal energy conversion flux is also determined by the density of the 
lower density stream. The conversion rate R decreases with increasing x, so the energy conversion 
is less efficient if the density contrast is large. In the limit of high x (x ~^ °°)? the collision of 
two streams reduces effectively to a problem where a solid body moving into a low density stream. 
Then, only one bow shock propagates into the low density medium and the postshock pressure is 
the ram pressure on the solid body. The contact surface moves with the velocity of the solid body. 

Another extreme is case the one where the colliding gas streams are identical (x = 1) and 
have opposite velocities, v 2 = —v\. Then the postshock regions have p 3 = p 4 = (4/3) / oif^ and 
v 3 = t> 4 = 0, and the shock speeds relative to the upstream flows are u\ = —u 2 = — (4/3)fi. The 
thermal energy conversion flux is F t h = (4/3) / oifj ! , and the conversion ratio is R = 4/3 which is 
the largest. The reason why R is greater than one can be understood as follows. All the kinetic 
energy of the gas that passes through the shock is converted into the thermal energy, since the 
postshock gas is at rest in the lab frame. But the shocks move outwards from the contact point, 
so R = |wi/fi| which is 4/3 for 7 = 5/3. 
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3. Three-dimensional Hydrodynamic Simulations 

We are interested in the situation where the gas streams move into an extremely low-density 
and low-pressure medium with a supersonic speed. For numerical simulations, however, we have 
adopted the flow parameters that can be handled comfortably with our code. In the collisions of 
two identical streams, we have set p = 1, p = 2.4 X 10~ 6 , and v = 0.2 for the streams, so the Mach 
number is v/c s = 100. Since the colliding gas is cold, the collision with a higher Mach number 
would represent a more realistic situtation. We have chosen the Much number 100, since it is high 
enough but still can be handled comfortably with our code. It is expected that the streams of stellar 
debris may have cross sections inversely proportional to the density, i.e., a ■ p = constant, while 
they have similar orbital speeds (Kochanek 1994; Monaghan & Lee 1994). Thus in the collisions 
of two different streams, we vary the density and the cross section according to a ■ p = constant. 
We use p 2 = 1 for the low density stream and the density contrast is x = /°i//°2 = (J 2/ (J i = 4, 9, 
or 16. The speed and pressure are the same, that is, v\ = v 2 = 0.2 and p 1 = p 2 = 2.4 X 10~ 6 . For 
the density contrasts considered here (x =4 to 16), we expect that the results would be somewhere 
between those of the identical streams (x = 1) and those of a high density contrast (x oo). 

The density and pressure of the ambient medium are set to be pf, = 10~ 3 pi, and pf, = 10 _3 pi. 
Since the pressure of the stream is higher than that of the ambient medium, a shock propagates 
into the ambient medium and a rarefaction wave into the stream. But the shock is weak with the 
Mach number w ps /c6 = 1.8, where c& is the sound speed of the ambient medium. If the ambient 
medium is a vacuum, the shock has the maximum speed u ps = 2cf ) /(j — 1). However, the maximum 
shock speed is still much smaller than the stream speed, so the shock propagating from the stream 
boundary is not dynamically important at all. 

For the boundary condition of the streams, we have used the inflow condition in which the flow 
variables of the streams entering the computational domain are continuous across the boundary 
(i.e., no gradient). Such inflow condition can be handled exactly, only when the inflow angle relative 
to the boundary is given by tan 9 m = ri\jn 2 where ri\ and n 2 are integers from to oo. Hence, in our 
simulations we consider only collisions with collision angles given by tan(#/2) = 1/3, 1/2, 1, 2, 3, 
and oo. The most likely collision angle between the outgoing and incoming gas streams expected 
around a black hole is greater than 90°; more restrictively between 130° and 140° as mentioned in 
§1. So we pay special attention for the case of 9 = 143° and 180°. 

We divide the simulations into two groups. Group A is for the simulations with two identical 
streams colliding with different collision angles. So Model "A36", for example, is for the collision 
of two identical streams with the collision angle 36°. Group B is for the simulations with two 
different streams colliding with either 143° or 180°. The stream parameters, such as the ratios of 
the densities and radii of the cross sectional areas of the two streams, are given in Table 1. 

The simulations have been done with a three-dimensional hydrodynamics code based on the 
Total Variation Diminishing (TVD) scheme (Harten 1983; Ryu et al. 1993). The TVD scheme is 
an explicit, second-order, Eulerian finite difference scheme. In the simulations with the TVD code, 
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shocks spread over 2 — 3 cells. In our simulations the computational box is a unit cube typically 
with 128 3 cells. Since the fast stream has a speed 0.2, it can cross the box in / = 5. The streams 
has been set to collide at / = and the temporal integration have been done up to i = 2.0. A 
typical run takes less than ~ 1000 CPU seconds on a Cray C90. 

4. Results and Implications 
4.1. Morphology 

Let us start by overviewing the morphological shape in the stream-stream collisions. To see 
some of the key features described in the following, refer to Figs. 2 to 5. Fig. 2 shows the pressure 
contour maps of a slice cut through the middle plane of the colliding streams for four models in 
Group A (A18, A36, A90, and A180). Fig. 3 shows the pressure contour maps for models of 9 = 143 
(i.e., A143, B143a, B143b and B143c). The density contour maps for the same models as shown 
in Fig. 3 are plotted in Fig. 4. Fig. 5 shows the velocity maps for models A143, B143c, A180, and 
B180. Morphologically the most noticeable feature is the oblique shocks formed around the contact 
surface. They can be identified by strong gradients in the pressure contour maps. 

For the collisions of two identical streams, one can define a symmetry plane which contains 
the initial contact surface. The angle between the oblique shocks and the symmetry plane depends 
upon the flow parameters and the collision angle. This angle is zero for the collision angle 9 = 180° 
and increases up to a maxim angle ~ 11° as 9 decreases to 90°. Then it decreases, as the collision 
angle further decreases below 90°. Once the angle between the oblique shocks and the symmetry 
plane is known, the shock jump condition can be found by applying the usual Hugoniot-Rankine 
relations. As the collision angle increases, the velocity component perpendicular to the symmetry 
plane increases so that the postshock gas pressure behind the oblique shock increases. 

The collision of two identical streams with the collision angle 9 = 180° (A180) is compared most 
closely to the one- dimensional collisions considered in the previous section. In the one- dimensional 
collisions, the shocks propagate back to the upstream flows and the shock speed relative to the 
upstream flows is 4/3f S f. But in the three-dimensional collisions the shock speed is ~ v st , since the 
shocks stand off near the symmetry plane instead of penetrating into the upstream flows. This is 
because the high pressure behind the shocks pushes the shock gas into the low pressure ambient 
medium in the direction perpendicular to the stream flows. The shocked gas is first accumulated 
into a dense region between the shocks, and then accelerated into the ambient medium. This 
expanding hot gas drives another shock in the ambient medium. This secondary shock is identified 
by strong pressure gradients in the background medium. The speed of this shock is similar to the 
stream speed itself, since the postshock pressure is similar to the ram pressure of the streams. The 
high pressure shocked region is roughly divided into a small volume of the densest gas surrounded 
by main shocks and a large volume of diffuse expanding flow encompassed by the secondary shock. 
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In the cases with the collision angle greater than 90° but less than 180° (e.g., A143), a forward 
shock and a reverse shock propagate into the two directions in the ambient medium. The forward 
shock moves faster than the reverse shock due to the positive bulk motion of the streams along 
the symmetry plane. The three-dimensional shape of the expanding high pressure region is then 
a torus that encompasses the contact surface, but the cross section of the torus is not symmetric. 
When the collisional angle is smaller than 90° (see Fig. 2), the reverse shock is absent. 

In the collisions of two different streams (B143abc and B180), the high density stream pen- 
etrates somewhat into the low density stream since the high density stream has a smaller cross 
section. However, the oblique shocks extend beyond the overlapping contact surface which is de- 
termined by the cross section of the high density stream. The flow outside the overlapping contact 
surface is prevented from streaming freely by the extended oblique shocks and the torus of the high 
pressure gas. Hence, contrary to a naive expectation, the geometric cross section of interacting 
streams seems to have a rather weak influence on the dynamics of the stream-stream collisions. 



The component of the orbital kinetic energy of the streams due to the motion perpendicular to 
the oblique shocks is converted to the thermal energy by heating the shocked gas. But the heated 
gas expands into the ambient medium and cools adiabatically. Through this adiabatic expansion, 
the thermal energy is converted back to the kinetic energy of the expanding gas encompassed by 
the secondary shocks, as described in the previous subsection. In order to quantify the efficiency of 
the conversion of the stream's orbital kinetic energy into the thermal energy, we define the energy 
conversion rate, R, as the ratio of the increase in the total thermal energy in the computational 
box to the kinetic energy added to the box through the boundary, i.e., 



where Ut t(t) is the total thermal energy at /. This ratio has a similar meaning to the conversion 
rate defined in Eq. (7). i?'s for various models are plotted against / in Fig. 6. The slow increase 
in R during the early phase for 6 less than 180° is due to the fact that the area of the contacting 
surface increases with time. The subsequent decrease is due to the following reason. The total 
thermal energy Ut t(t) increases initially with /, but then stays more or less constant after the 
contact surface is fully developed. But the kinetic energy in the denominator of Eq. (9) continues 
to increases with time, so the ratio R decreases at late epochs. In the end of our simulations (/ = 2), 
only a small fraction of the stream's orbital kinetic energy (^30%) was converted to thermal energy 
and most of it went to the expanding kinetic energy of the hot torus gas. 

Since the initial orbital energy of the streams is converted eventually to the expanding energy 
of the shocked gas, the nature of the debris orbit would change drastically after the stream crossing. 
The amount of energy transferred to the expanding motion can be estimated from the maximum 



4.2. Energy Conversion 



R = 




(9) 



(l/2)p 1 vl(a 1 v 1 t) + (l/2)p 2 vl(a 2 v 2 t) j 
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value of the energy conversion rate. If the maximum R is close to 1, almost all the stream's orbital 
energy will be eventually converted into the expanding energy. In the cases of the collisions of two 
identical streams, the maximum R is generally close to 1. The maximum R becomes smaller as 
the cross sectional area (and thus density) contrast grows, but the dependence is rather weak. For 
example, the maximum i?'s of the different stream cases (for x = 4 — 16) are about 30 — 35% of 
that of the identical stream case (see Fig. 6c). Thus the conversion of the orbital energy to the 
expanding energy is still significant unless the density contrast is very large (i.e., x ^ 10). This 
is mainly due to the fact that the oblique shocks affect the area beyond the overlapping contact 
surface, as discussed in the previous subsection. Similar level of decrease in the energy conversion 
with the density contrast is seen from Fig. 6d (6 = 180°) where head-on collision cases are shown. 



The simulations of colliding two streams would comprise of only a small part of the evolution 
of stellar debris around a black hole, but we may still extract some important physics. In this 
subsection, we attempt to predict how the collision of debris streams affects their evolution and 
how the circularization of the debris orbit proceeds, by relating the results of the present simulations 
to those of stellar disruption simulations described in Monaghan & Lee (1994). 

In Fig. 7, we show the typical geometry of gas streams resulted from an SPH simulation for the 
disruption of a star near a black hole (Monaghan & Lee 1994). The arrows indicate the direction 
and the size of the velocity of SPH particles. The magnitude of the velocity of the streams is close 
to escape velocity, because the eccentricity is almost 1. Therefore, at a distance r from the black 
hole it can be approximated as 



where M is the mass of the black hole. The angle of the stream crossing is approximately given as 



where r p is the pericentral distance of the orbit of the disrupted star and Rs is the Schwarzschild 
radius of the black hole. Here we have assumed that the velocity vector at the crossing is nearly 
radial. The relation given by Eq. (10) can be confirmed from Fig. 8, where we have shown the 
velocity of SPH particles as a function of distance from the black hole. Note that the length unit 
shown in Fig. 7 and 8 is Rs/2. The circles denote the outgoing particles while the crosses represent 
the incoming particles. There is no distinction in orbital velocities for the incoming and outgoing 
particles. Since the location of stream crossing is about r cr = 50Rs (or 100 in present length unit) 
from the black hole, it can be seen that the velocities are closely approximated by Eq. (10). With 
a typical value of r p ~ r t id ~ 5Rs, we obtain 9 ~ 130°, where r t id is the tidal disruption distance. 



4.3. Implications on the Long Term Evolution of Stellar Debris 




(10) 




(11) 
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We can estimate the average velocity of the shocked streams along the tangential direction as 
v tan ~ v st \/^ cos(#/2) for the collision of two identical streams, if we assume that the collision 
is sticky (i.e., inelastic) and the gas is incompressible. It is equal to the circular velocity of 
v c i r = \/GMjr for 9 ~ 151°. Thus for the typical angle of 9 ~ 130°, the tangential velocity 
of the sticky, incompressible fluid will be close to the value for a circular orbit. In such case the 
orbital energy would be removed effectively and the shocked gas would move nearly circular orbit 
at the orbital radius r cr . In reality, however, the gas is not sticky nor incompressible, so the shocked 
gas goes through an adiabatic expansion. 

The geometry of the crossing debris streams would resemble that somewhere between A143 
and B143c (see Figs. 3 and 4). In Fig. 9 we show the accumulated frequency distribution of the 
kinetic energy per unit mass (Ek) for all the gas within our simulation box at the end of our runs of 
models A143a and B143c. The kinetic energy of the initial streams (Ek )S ) 1S indicated as a vertical 
line at Ek/ ' Ek )S = 1. The presence of gas at both lower and higher sides of the initial Ek is due 
to the shock. Bulk of the shocked gas have energy smaller than the initial stream's kinetic energy. 
Thus the stream crossing will make the shocked gas more tightly bound to the black hole. A small 
fraction of the shocked gas has kinetic energy greater than that of the initial streams, as evident 
from the high velocity tail in Fig. 9. The amount of gas having Ek > Ek )S is about 1/4 of that with 
Ek < Ek tS . However, still majority of the high velocity gas lies very close to Ek s , so the fraction of 
gas that will be unbound from the black hole is likely to be very small (i.e., < 10%). 

The difference in the circularization efficiency between the identical and different stream 
collisions can be noticed also from Fig. 9. About half of the gas with Ek < Ek )S have Ek < 0.5Ek tS 
for the model A143a, while the fraction becomes less than a third for the model B143b. Ek = 0.5Ek tS 
represents approximately the energy of the gas having a circular orbit after the shock. Thus the 
particles with Ek < 0.5Ek tS have a apocenter smaller than r cr , while those with E > 0.5Ek tS have 
a apocenter greater than r cr . 

Once the stream crossing begins, the incoming stream no more exists beyond the crossing 
point. The tail of the incoming stream moves around the black hole and becomes the endpoint 
of the outgoing stream. After the tail of the outgoing stream is destroyed by the collisions, the 
incoming stream will flow freely again. Thus the duration of a stream-stream collision until the 
exhaustion of one segment of the outgoing stream is about the time required for a particle to have a 
round trip from the location of the stream crossing through the apocenter from the black hole. For 
a particle with a nearly parabolic orbit, the duration is twice the free fall time. Thus the collision 
time scale would be t co \ ~ 2^/ (3ir 2 r 3 ) / (8G M) ~ 50 hours for our choice of typical parameters (i.e., 
M = 1O 7 M , r cr « 50^5, 

and r, p ~ r t id ~ 1 AU). The quiescent phase would also have the similar 
duration. Therefore, the stream collision could be a periodic phenomena with period about 100 
hours. A tidal disruption with smaller r p would produce shorter period, because r cr is smaller. 
Obviously, such a periodicity is less pronounced if the contrast in cross sectional areas becomes 
larger. 
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5. Conclusion 

In this paper, we have investigated the supersonic collisions between two gas streams with 
various configurations and collision geometry. In the collisions of one- dimensional plane-parallel 
streams, the kinetic energy is converted effectively into thermal energy via the shocks which 
propagate into the streams. The conversion is most effective for the collisions of two identical 
streams. It becomes least effective if the two streams have an infinite density contrast. For the 
intermediate density contrast, the thermal energy conversion is proportional to the ram pressure of 
the low density stream but nearly independent of that of the high density stream. If the shocked 
gas can cool radiatively in a time scale shorter than the dynamical time, the kinetic energy of the 
streams would be dissipated more effectively. 

The collisions of three-dimensional streams are much more complicated. We have more degrees 
of freedom in setting flow parameters such as the ratio of cross sections and collision angle. The 
key features of the three-dimensional simulations which contrast with those of the one- dimensional 
collisions include: 1) The shocks generated by the collisions stand near the stand-off point instead 
of penetrating into the colliding streams. 3) Most of the bulk kinetic energy of the streams turns 
into the kinetic energy of the expanding hot gas instead of the thermal energy. 2) The hot gas 
behind the shock expands uniformly in the directions perpendicular to the colliding streams. 

In the collision of two identical streams, the fraction of the orbital kinetic energy converted to 
the expanding kinetic energy is high, so the circularization is effective. But realistic streams at the 
crossing point are likely to have a density contrast as large as ~ 4 — 10 and the ratio of cross sections 
is expected to vary as the reciprocal of the density contrast. In such cases, the conversion of the 
orbital kinetic energy should be somewhat smaller. If we quantify the circularization efficiency by 
the fraction of the orbital kinetic energy transferred to the expanding kinetic energy, the efficiency 
would be nearly 50% for B143a (x = 4) and B143b (x = 9) and 30% for B143c (x = 16), while 
nearly 100% for A143 (x = 1)- The reason that the circularization efficiency depends only weakly 
on the ratio of cross sections is that the shock heated gas forms a thick torus of high pressure gas 
which keeps the stream of larger cross section from streaming freely. Thus stream-stream collision 
could still be an effective circularization mechanism, even if the ratios of densities and cross sections 
of two streams are significant (x^lO). 

The detailed evolution of the expanding gas under the gravity of the black hole is still com- 
plicated. Some would escape from the black hole gravity, taking energy away with it and leaving 
the debris more tightly bounded to the black hole. Therefore, the collision could be considered as 
an effective way to dissipate the orbital kinetic energy of the streams, even in the absence of the 
radiative cooling of the shocked gas. 

The stream crossing would last for a while (about 50 hours for typical parameters, i.e., 
M = 1O 7 M0, r cr ~ 50i?5, and r p ~ rm ~ 1 AU) with a similar period of quiescent phase followed. 
Such a pattern is likely to repeat until the majority of stellar debris returns to the apocenter. Thus 
a periodic occurrence of flare could be observed if the stream crossing can produce any observable 
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phenomena. 

So far, we have considered an adiabatic gas with the adiabatic index 7 = 5/3 and ignored 
any radiative processes. However, in the collision of two gas streams, their kinetic energy could 
go mostly into radiation pressure rather than gas pressure depending on the parameters. Then 
7 = 4/3 would be more appropriate than 7 = 5/3. But we believe that this is unlikely to make 
much of a difference to the outcome of our calculations. Also, it is possible that radiative processes 
might be important. For instance, if the torus of the shocked gas is optically thin, radiative losses 
could be efficient and remove internal energy almost as fast as it is created in the shocks. If it is 
optically thick, then radiative viscosity might play an important role. But since we have limited 
ourselves to the dynamics of an adiabatic gas in this paper, we leave these to future work. 
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Table 1. Model Parameters 



Model 





P1/P2 


ri/r 2 


A18 


18 


1 


7/7 


A26 


26 


1 


7/7 


A36 


36 


1 


7/7 


A45 


45 


1 


7/7 


A52 


52 


1 


7/7 


A63 


63 


1 


7/7 


A71 


71 


1 


7/7 


A90 


90 


1 


7/7 


A126 


126 


1 


7/7 


A143 


143 


1 


7/7 


A180 


180 


1 


7/7 


B143a 


143 


4 


5/10 


B143b 


143 


9 


3/9 


B143c 


143 


16 


3/12 


B180 


180 


16 


3/12 
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Fig. 1. — Flow structure after the collision of one- dimensional plane-parallel gas streams. The flow 
parameters pi,pi,v\, and P2,P2, v 2 represent the initial condition left and right to the discontinuity. 
Two shocks propagate away from the discontinuity. 

Fig. 2. — Pressure contour maps of a horizontal cut in models A18 (top left), A36 (top right), A90 
(bottom left), and A180 (bottom right). 

Fig. 3. — Pressure contour maps of a horizontal cut in models A143 (top left), B143a (top right), 
B143b (bottom left), and B143c (bottom right). 

Fig. 4. — Density contour maps of a horizontal cut in models A143 (top left), B143a (top right), 
B143b (bottom left), and B143c (bottom right). 

Fig. 5. — Velocity field maps of a horizontal cut in models A143 (top left), B143c (top right), A180 
(bottom left), and B180 (bottom right). 

Fig. 6. — Energy conversion rate R. Top two panels are for the models with two identical streams. 
R increases with the collision angle. Bottom left panel is for the models with the collision angle 
143°. The solid line is for A143 (i.e., identical streams), the dotted line for B143a (x = 4) the short 
dashed line for B143b (x = 9), and the dot-dashed line for B143c (x = 16). Bottom right panel is 
for the models with the collision angle 180°. The solid line is for A180 (i.e., identical streams), the 
short dashed line for B180 (x = 16). 

Fig. 7. — A typical geometry of stream crossing resulting from an SPH simulation. The arrows 
indicate the locations and velocities of SPH particles. The length unit in this figure is Rs/2. 

Fig. 8. — The velocity of SPH particles as a function of distance from the black hole for the case 
shown in Fig. 7. The straight line indicates the v — r relationship for a parabolic orbit. 

Fig. 9. — Accumulated frequency distribution of the kinetic energy per unit mass (Ek) for the 
entire gas in our simulation box at the end of simulations for models A143 and B143c. The kinetic 
energy of the incoming stream is indicated as a vertical line at E^/ Ek tS = l. Majority of shocked 
gas has lost its kinetic energy while a small fraction has gained it. 
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